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Abstract 

Since the calculation of BR(B~ — * rj'K") in the framework of 
QCD improved factorization method, developed by Beneke et al, leads 
to numerical values much below the experimental data, we include two 
different contributions, in an alternative way. First, we find out that 
the spectator hard-scattering mechanism increases the BR value with 
almost 50%, but the predictions depend on the combined singularities 
in the amplitude convolution. Secondly, by adding SUSY contribu- 
tions to the Wilson coefficients, we come to a BR depending on three 
parameters, whose values are constrained by the experimental data. 



1 Introduction 

As a first evidence of a strong penguin, the B~ — > rfK~ decay has be- 
come of a real interest after CLEO announced its large numerical value 
BR(B~ -> rfK~) = (6.5t};f ±0.9) x 10~ 5 [1], which could not be explained by 
the existent theoretical models. As improved measurements followed, provid- 
ing even larger values, (80^°±7) x lO" 6 (CLEO [2]), (76.9±3.5±4.4) x lO" 6 
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(BaBar [3]) and (79±u ± 9) x 1(T 6 (Belle [4]), the inclusion of new contri- 
butions for accommodating these data has quickly become a real theoretical 
challenge. In this respect, perturbative QCD mechanisms [5], with different 
rj'g*g* vertex function [5, 6], have been considered as main candidates for sig- 
nificantly increasing the BR(B~ — > rj'K~) value. On the other hand, while 
searching for physics beyond the Standard Model (SM), supersymmetry has 
been employed in processes like B — > J/ipK* [7], B — > <f)K [8], B — > nK [9, 
10], B — > X s 7 [11], and deviations from the SM predictions for the values of 
branching ratios and CP asymmetries have been the main targets. 

The present paper is organized as follows: in Section 2, we compute 
the BR(B~ — > r)'K~) in the improved factorization approach developed by 
Beneke et al. [12]. Since we get a BR much below the experimental values, 
we incorporate two alternative contributions. The first one, presented in 
Section 3, comes from the so-called spectator hard scattering mechanism. 
Following a similar approach as in [13], we give a detailed calculation of the 
gluonic transition form factor which plays an important role in the evaluation 
of this contribution. Although it has been concluded that this mechanism 
could provide large BR values [13], we show that the presence of combined 
singularities in the amplitude convolution is a source of large uncertainties. 
In Section 4, we employ a supersymmetric approach and include exchanges 
of gluino and squark with left-right squark mixing. Working in the mass 
insertion approximation [14], the values of the Wilson coefficients c$ g and c 77 
can be significantly increased, by adding the SUSY contributions, and this 
has a strong numerical impact in the branching ratio estimation. Finally, one 
may use the experimental data to impose constraints on the flavor changing 
SUSY parameter 5 b L s R . 
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2 Improved QCD Factorisation 



The relevant decay amplitude for B — > rj'K , in the improved QCD factor- 
ization approach [12], is given by [5, 15] 

A{B-^n'K~) = -i^(m|-mJ)Fo BW (<)/^[K 6 K>iW 

+ V pb V; s (a p 4 (X) + a p 10 (X) + rf (a£(X) + af (X)))} 

- i^=(m 2 B - m 2 K )F*- K {ml)f% [V ub V: s a 2 (Y) 
+V pb V p * s [(a 3 (Y) - a 5 (Y)) (2 + a) 



<(Y)--aUY)+r' x U{Y) - -a*(Y) 



+ 



+ -(a 9 (Y)-a 7 (Y))(l-a) 



a 



(1) 



where X = rj'K and Y = Krj', p is summed over u and c, r' x = 2m^,/(m b 
m s )(2m s ), r* = 2m 2 K /m b (m u + m s ), a = f°,/f%, and [12] 
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where Cf = (N% — 1)/2N C and N c = 3 is the number of colors. The vertex, 
the hard gluon exchange with the spectator and the penguin contributions, 

at (j, — rrib, are: 
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where x = 1 — x, y = 1 — y and the parameter 2ji M /m h coincides with r x . 
The functions g(x), Gm{x) and G\m(x) are given by 
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Gm(s) = [ dxG(s — ie,x)4>M(x), 
Jo 

G M (s) = / dxG(s - ie,x)(f) p M (x), (4) 
Jo 

where Sj = mf/ml are the mass ratios for the quarks involved in the penguin 
diagrams, namely s u = Sd = s s = and s c = (1.3/4.2) 2 . 

As it can be noticed, except for the hard contribution where the wave 
functions for both Mi and M2 are involved, the coefficients a>i are different 
for the X and Y final states, since they depend on the twist-2 and twist-3 
wave functions of the M2 meson. Thus, the twist-2 distribution amplitude 
4>k(x) has the following expansion in Gegenbauer polynomials [12, 16] 

<f) K (x) = 6x{l -x)[l + af C[ m (2x - 1) + af Cl' 2 (2x - 1) + ...], (5) 

with C? /2 (u) = 3m, Cl ,2 {u) = (3/2) (5m 2 - 1), af = 0.3 ± 0.3, and of = 
0.1 ± 0.3. The corresponding twist-3 amplitude, <p p K) is 1. 

The physical states 77 and 77' are mixtures of SU(3)-singlet and octet com- 
ponents i] and i] 8 and therefore the corresponding decay constants, in the 
two-angle mixing formalism, are given by 

^ = vl sin 08 + {k cos 60 ' 

f°, = - 2 A sin fl 8 + _| cos £ 0) (6 ) 

with ^ 8 = -22.2°, 9 = -9.1°, / 8 = 168 MeV, and f = 157 MeV [17]. These 
lead to /"/ = 63.5 MeV, f*, = 141 MeV and to the relevant form factor for 
the B — > r]' transition 

^r^w ^/sin^a cos#n\ / \ 

F o = + -^J = 0.137 (7) 

Even though the r( flavor singlet meson has a gluonic content which could 
bring a contribution to the wave function, this is supposed to be small [18] 
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and therefore we employ, in the calculation of V v >, P^ 2 and P^' 2 in Oj(Y), 
only the leading twist-2 distribution amplitude 

4>r)' = Gxx . (8) 

Also, since the twist-3 quark- ant iquark distribution amplitude do not con- 
tribute, due to the chirality conservation, the penguin parts in a p & {Y) and 
Og(y) are missing. As for the B meson wave function, we shall work with 
a strongly peaked one, around z = \ B /m B ss 0.066 ± 0.029, for \ B = 
0.35 ±0.15 GeV. 

Putting everything together, we get, within the SM improved factoriza- 
tion approach [12], the numerical value BR S m(B — > rj'K) = 3.65 x 10" 5 , 
which although is in accordance with other theoretical estimations [5, 15, 
17], yet lay below the experimental data [1-4]. Hence, in spite of the "conser- 
vative" prediction that the conventional mechanism should be the dominant 
one, it has been getting clear that new contributions are needed in order to 
account for the existent data. 



3 Spectator Hard-Scattering Mechanism 

It has been considered that the spectator hard-scattering mechanism (SHSM), 
depicted in Figure 1, is a reliable framework for this process, which signif- 
icantly increases the value of BR(B — > rj'K) [5, 13]. Following this idea, 
let us write down the corresponding di-gluon exchange amplitude for the b 
quark decaying into an s quark and a hard gluon 
/b Ik f 

A hs = -iC F gl B K J dzdy(j) B (z)(f) K (y) 



x Tr [ lb P k V,{P B + m B ) l5lu ] (Q?, Qj, m%) (9) 



in terms of the effective b — > sg vertex [19] 
T a = Gf g. 
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will play an important role in the evaluation of the amplitude A^. Perform- 
ing the calculations in (9), we come to the following expression of the hard 
scattering amplitude: 

A hs = -2t^v; s v pb ^f B f K (2f^ + f^jy z( t )B ^)l 1 dyMy) 
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and neglecting, for the moment, both m 2 , and m 2 K , the amplitude (14) be- 
comes 
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where, for the dominant contribution coming from the insertion of the O"' 
and the magnetic-penguin 0$ g operators, one has [13] 
"2 
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In what it concerns the F(y, a) function, which is an essential input in the 
calculations, it can be first written as 

6xx (Ql + Q\ — 2xxm 2 l ,) 
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where we have introduced the notation p 2 = Ql + Q 2 — m 2 , . The logarithmic 
nature of the F(y,a) function makes it very sensitive to the values of Ql, 
Q 2 ,, m 2 ,. We recommend [6] for a detailed discussion of the f]'g*g* vertex 
in the case of arbitrary gluon virtualities in the time-like, Ql > 0, Q\ > 0, 
p 4 - AQlQl > 0, and space-like, Ql < 0, Q\ < 0, p A - AQlQ 2 < 0, regions. 
Now, using 
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where we have neglected m 2 K , the dominant term in (19) is: 
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On the other hand, by comparing the expressions in (20), it clearly results 
that we are in the limit where \Ql\ ^> \Q\\- So, the function F(y,a) can be 
computed in this approximation and it simply yields 



F{y,a) = T 

m 2 , 



y- z 

a 2 



y) log 



l - 



(22) 



As it can be seen from (20), the term (y — z )/a 2 + y = Q\/m 2 rj , takes a 
whole range of values, from —0.87 to 26.5, as Q\ goes from the space-like 
to the time-like regions. Consequently, a logarithmic singularity develops 
as y — > zq/(1 — a 2 ), i.e. for Q\ — > m 2 ,. Inspecting (16), we also notice 
the pole at y = in the contribution. In addition, while G[s p , (1 — 
z )(y — zq)\ is divergence free for all s > 0, the G[0, (1 — z )(y — z )] gets 
a logarithimic singularity at y = z . Hence, in the course of numerically 
evaluating the scattering contribution, one must be careful about dealing 
with these combined singularities in the convolution (16). 

As in the case of other hard-scattering theoretical estimations [5, 13], the 
amplitude of this contribution contains, as main uncertainty, the peaking po- 
sition, Zq, in the B meson distribution function and accordingly, the branch- 
ing ratio is extremely sensitive to it. For z e [0.063, 0.068] and the average 
value a s (Q\) = 0.28, the total branching ratio, including besides the im- 
proved factorization approach, the spectator hard-scattering mechanism with 
the vertex function (22), is in the range from BR(B — > rj'K) = 6.58 x 10~ 5 , 
for z = 0.063, to BR(B -> -q'K) = 5.8 x 10~ 5 , for z = 0.068. 

Comparing these results with the experimental data [1 — 3], we notice 
that they are still below the lowest limit. An alternative way which increases 
the BR and avoids the uncertainties coming from the combined singularities 
in the convolution (16), would presumably look more reliable. 

4 SUSY Gluonic Dipole Contribution 

Employing the Minimal Supersymmetric Standard Model (MSSM), we shall 
add to the effective SM Hamiltonian (1), the SUSY contribution 

H SUSY = - i %iV ub V: s + V cb V* s ) (c!^0 8g + c^0 77 ) , (23) 
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expressed in terms of the usual gluon and photon operators: 

°s 9 = ^m b sa^(l + l5 )G^b, 

7l = -^m 6 ^(l + 75 )F^6, (24) 



and of the Wilson coefficients [10, 20] 
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where 



G (x) = X [22 - 20x - 2x 2 + 16x \n(x) - x 2 ln(x) + 9 \n(x) 

4:X r 1 

F (x) = - -r: 1 + 4a; - 5x 2 + Ax ln(x) + 2x 2 \n(x) (26) 

In the above expressions, x = m 2 l /m'~, with trig being the gluino mass and trig 
an average squark mass, while the factor 5 bs = A 6s /m~, where A bs are the 
off-diagonal terms in the sfermion mass matrices, comes from the expansion 
of the squark propagator in terms of 5, for A <C m~. In principle, the 
dimensionless quantities S bs , measuring the size of flavor changing interaction 
for the sb mixing, are present in all the SUSY corrections to the Wilson 
coefficients in (1) and they are of four types, depending on the L or R helicity 
of the fermionic partners. In the followings, we focus on the 8 b £ R insertions 
because only the SUSY Wilson coefficients (25), being proportional to the 
large factor m g /m b , are going to make an important contribution, even for 
small values of 5. 

In (3), we replace the Wilson coefficients c%J and cfif , by the total quan- 
tities 

C [*> 5 ] = 4 f + 4e SY (™ b ) , c% al [x, 5} = c% + c s ^ S¥ (m b ) , (27) 
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where c SUSY (m b ) have been evolved from M SUSY = m g down to the \i = m h 
scale, using the relations [10, 19] 

c S ^ SY {m b ) = ^^(m^ + ^^-^cf^K), (28) 



with 



7/ = (a s ( m 9)/«s( m <)) 2/21 (« s (m t )/a s (m fe )) 2/23 (29) 



We choose for rriq the value rriq = 500 GeV and write m g as = ^fx rriq 
and 5^, = pe lip . As the total branching ratio can be expressed in terms of 
three free parameters: x, p, <p, one is able to plot the BR total , in units of 10 -5 , 
as a function of (p, ip), for different values of x. By inspecting the 3D plots 
displayed in Figure 2, for x = 0.3 (the upper) and x — 1 (the lower surface), 
we notice that the SUSY contributions (25) to the Wilson coefficients have 
significantly increased the SM value, BR SM = 3.65 x 10~ 5 , represented by the 
horizontal plane. Using the experimental data, one is able now to determine 
the 8 h L R complex values, for each x. 

Let us take, for example, x — 1, pointing out that the same discussion 
can be performed for any rr-value. For p = 0.005, the BR total is increasing 
from 5.1 x 10~ 5 , for ip pa ±7r/3, to the maximum value BR total = 6.24 x 10~ 5 , 
for tp — 0. As p goes to bigger values, we find a better agreement with the 
large experimental data. For p = 0.01, the data can be accommodated for 
ip pa — 7r/4, while, for p = 0.02, one has to impose <p pa — Sn/15. 

5 Concluding Remarks 

At first, we have analyzed the B~ — > rj'K decay and computed its branch- 
ing ratio using the improved factorization method developed by Beneke al. 
[12]. Since the obtained result, BR S m = 3.65 x 10~ 5 , is much below the 
experimental data, [1-4], we have have added new contributions. 
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In this respect, the so-called spectator hard-scattering mechanism, which 
is depicted in Figure 1, has allowed us to compute the amplitude in terms of 
the effective b — > sg vertex and the transition form factor (11) which contains 
the quark contribution to the r)'g*g*, (21), as an essential input. The total 
BR has, as a main uncertainty, the peaking position in the B meson wave 
function, z = A B /m B , with \ B = 0.35 ± 0.15 GeV. Even the results are 
closer to the experimental data, we point out the combined singularities in 
the amplitude convolution (16) which must be treated carefully. 

Secondly, we extend the SM to the MSSM and add SUSY contributions 
to the Wilson coefficients cj g and c 7 lf . The total BR is expressed in terms 
of the parameters x = m^/m^, and 5 b ^ R = pe %Lp whose contribution turns out 
to be important, even for very small values of p. Finally, by inspecting the 
3D-graphics (see Figure 2), representing the BR total for x = 0.3 (the upper 
surface) and x = 1 (the lower surface), one is able to find numerical values 
for p and (p that can account for the experimental data or other theoretical 
predictions [21]. 
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Figure Captions 



Fig.l. Feynman diagrams of the hard scattering mechanism for B — > rj'K . 
The gluons are represented by the dashed lines. 

Fig.2. Total branching ratios (SM+SUSY) for B~ -> rfK~, in units of KT 5 , 
as functions of (p, ip), for x = 0.3 (the upper plot) and x — 1 (the lower plot), 
compared to the SM estimation represented by the horizontal plane. 
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